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ABSTRACT 

The aim of this paper is the verification of the attenuation models: 
Atkinson&Boore 2008, Bindi et al. 2013, and Akkar & Bommer 2013, for cortical 
earthquakes registered in Quito city,through the RED-RENAC-IGEPN network. The 
systematization of the procedure described in the movement prediction equations, 
is done in Matlab programming language with the GUI Guide tool. For the handling 
and visualizaron of the results, the signáis are transformed into parameters for use 
engineering, like spectral accelerations, using the Newmark’s method. The Fourier 
transform is used to verify the signal frequency content, Scherbaum (2010) criteria 
are used for selection and validation described GMPEs. 

Keywords: Accelerograms, Response Spectra, Design Spectra, Attenuation 
Equations 

ECUACIONES DE ATENUACIÓN COMPATIBLES A LA 
TECTÓNICA DE QUITO GMPEs. 

RESUMEN 

El objetivo del presente trabajo es la verificación de los modelos de 
atenuación: Atkinson&Boore2008, Bindi et al. 2013, y Akkar&Bommer 2013, para 
los sismos corticales registrados en la ciudad de Quito, mediante la red RED 
RENAC-IGEPN. La sistematización del procedimiento descrito en las ecuaciones de 
predicción del movimiento, se realizó con el lenguaje de programación matlab, con 
la interfaz gráfica Guide, para el ingreso y visualización de sus resultados, luego las 
señales se transformaron a parámetros de uso ingenieril, como las aceleraciones 
espectrales, utilizando el algoritmo de Newmark. 

Para verificar el contenido de frecuencias de señales se utilizó la 
transformada rápida de Fourier. En tanto que para la validación de los modelos de 
predicción se tomó los criterios de Scherbaum (2010), mediante el índice log- 
Likelihood. 

Palabras Claves: Riesgo Sísmico, Acelerogramas, Espectros de respuesta, 
Espectros de diseño, Ecuaciones de atenuación. 
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1 INTRODUCCIÓN 

En el Ecuador el origen de los sismos provienen de dos fuentes: La primera 
por subducción de la placa de Nazca, que induce a la placa Sudamericana y la 
segunda, debido a la presencia de fallas locales corticales, denominadas sistema 
de fallas de Quito (QFS). 

En el caso de la ciudad de Quito, la desagregación de la amenaza muestra 
que los sismos de mayor intensidad provienen de las fallas corticales (Alvarado et 
al. 2014), por esta razón, para aportar al mejoramiento de los cálculos probabilistas 
de la amenaza, tanto en los estudios de la normas sísmicas, como en proyectos 
ingenieriles, es necesario acudir a la información de sismos de baja intensidad 
registrados por la Red Nacional de Acelerógrafos (RENAC). 

La Norma Ecuatoriana de Construcción (NEC-15), en su primera emisión 
año 2011, Apéndice 10, sección 10.1. “Fundamentos del Estudio de Peligro 
Sísmico”, afirma: “Las ecuaciones de predicción (GMPEs), se fundamentaron en 
modelos de predicción para fuentes intraplaca como: Boore-Atkinson (2008) y 
Akkar-Bommer (2010)”. El presente trabajo pretende verificarla aplicación de los 
modelos Bindi et al. (2013) y Akkar&Bommer (2013), según clasificación mostrada 
en la figura 1. 



Figura 1 . Modelos de predicción de GMPE’s, válida para la ciudad de Quito 


2 SISTEMA DE FALLAS CORTICALES EN LA CIUDAD DE 
QUITO (QFS) 

Los sismos en la ciudad ocurren principalmente por el sistema de fallas de 
Quito (QFS), conocidas como: “Puengasí”, “llumbisí-La Bota” y “Carcelén-EI Inca”, 
hacia el norte el sector de Guayllabamba se tiene los segmentos de falla “Bellavista- 
Catequilla”, “Tangahuilla, R.Uravia F”. y “RCoyango F (GFS)”. 
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La formación de estas fallas inversas se debe al movimiento del bloque 
formado por los valles de “Tumbaco” y “los Chillos”, los cuales forman un ángulo 
agudo, con el bloque de Quito, que al desplazarse generan esfuerzos cortantes. 
(Alvarado 2014). 

La formación este plano de falla se debe al movimiento relativo entre un 
bloque 1 (Hanging wall), respecto al bloque 2 (Foot wall), generando esfuerzos 
cortantes en este plano, como se muestra en la figura 2. 



Figura 2. Falla inversa en la ciudad de Quito. 


Dónde: el bloquel, es el que se desplaza (“Los Valles”-Foot wall), en tanto 
que, el bloque2, permanece fijo (“Quito”-Hanging wall) 

En la tabla 1, se muestra los sismos ocurridos en la ciudad de Quito, entre 
los años 2014 al 2016, y registrados por la Red Nacional de Acelerógrafos (RENAC), 
perteneciente al Instituto Geofísico de la Escuela Politécnica Nacional (IGEPN), la 
profundidad y el tipo de fallas se toma de los reportes emitidos por esta entidad. 
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TABLA 1 Sismos locales registrados por red RENAC-IGEPN 


Sismo 

Fecha 

Mw 

Depth 

(Km) 

dt 

Tipo de 
falla 

Calderón 

12-08-2014 

5.1 

5.3 

0.01 

Inversa 

Puembo 

8-08-2016 

4.6 

7.3 

0.01 

Inversa 

Puembo-Tab. 

4-09-2016 

4.7 

4.7 

0.01 

Inversa 

Conocoto 

17-02-2011 

4.2 

12 

0.01 

Inversa 

Guayllabamba 

29-10-2011 

4.1 

3 

0.01 

Inversa 

Calderón 2 

16-08-2014 

4.7 

7.7 

0.01 

Inversa 


3 MODELO DE WELLS Y COPPERSMITH 1994 


Mediante las ecuaciones de Wells y Coppersmith 1994, se calcula la 
longitud de la falla sub superficial (RLD) y el ancho de falla (RW). 



Figura 3. Representación del cálculo de las distancias (Conocoto 17/ 02/ 2011.) 
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En la tabla 2, se muestra las distancias, desde las estaciones instaladas, 
hacia el sitio de la falla con el modelo de Wells y Coppersmith 1994, a escala real 
mostrado en el esquema tipo de la figura 3, para los demás sismos descritos en la 
tabla 1, se procede de manera similar. 


TABLA 2 Distancias epicentrales y las estaciones 


Estación 

Repi. 

Rjb. 

Rhyp. 

PGA-NS 

PGA -EW 


(Km) 

(Km) 

(Km) 

(g) 

(g) 

LILI 

5.14 

5.38 

13.49 

0.014 

0.017 

AEPN 

4.6 

4.43 

13.2 

0.016 

0.02 

PRAM 

12.04 

11.65 

17.3 

0.043 

0.036 


Dónde: 

Repi.= Distancia epicentral (Km) 

Rjb.= Distancia Joyner Boore (Km) 

Rhyp.= Distancia hipocentral (Km) 

PGA= Aceleración pico del terreno de la estación (g) 

4 FORMA TÍPICA DE LAS ECUACIONES DE ATENUACIÓN 
GMPES 

La forma típica de la ecuación de predicción del movimiento fuerte es: 

Y = C l+ C 2 M + C 3 M c *+C 5 Ln[R + C 6 exp (C 7 M )] + C,R + f fuente + f sitio + <r (1) 
Donde cada término de la ecuación representa: 

Y = Valor pico del movimiento fuerte 

C 2 M + C 3 M Cí = Factor referente a la magnitud del sismo. 

C 5 ln [/?...] + C S R =Factor referente a la distancia. 

f 'fuente + fsitio = lnf| uencia del tipo de fuente y sitio. 

Cl,C2,...,Cn = Constantes de regresión 
<J = Desviación estándar. 

5 MODELO DE ATKINSON & BOORE 2008 (AB08) 

Este modelo predice el movimiento fuerte del terreno (GMPE’s) en función 
de la magnitud de momento, la distancia horizontal a la proyección del plano de falla 
denominado distancia Joyner-Boore, rango de la velocidad de onda de corte a los 
30m de profundidad, y el tipo de fallamiento. La ecuación determina la aceleración 
(PGA) y velocidad (PGV) pico del terreno, así como el espectro de respuesta (PSA, 
al 5% amortiguamiento pseudo-aceleración), para periodos comprendidos entre 
O.OIs y lOs. 

ln(Y) = F m (M) + F d {R ¡b ,M) + F S (V S30 ,R JB ,M) + ea T 


( 2 ) 
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5.1. Características del modelo (AB08) 

Este modelo es válido para el rango de magnitud de Mw: 5-8 grados, la 
distancia entre 0-400km, distancia Rjb, debe ser menor a 200km, Vs30 (180-1300) 
m/s, periodo 0.001-lOseg. 

Las componentes de sismo se calcula con la media geométrica los registros 
seleccionados, además, no es recomendable para suelos duros, con velocidad de 
onda de corte Vs30>1500m/s. 


ATKINSONKBOORE 2008 vs SEÑAL SELECCIONADA 
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Figura 4. Interface programa GMPFts-EPNvl, modelo de Atkinson & Boore 2008. 


En la figura 4, se muestra el espectro de respuesta, con la validación del 
programa GMPE-EPN.vl, de acuerdo a los siguientes datos: Magnitudes Mw=5, 6 
y 7, distancia Ftjb=200km, velocidad de onda de corte Vs3o=760m/s, y tipo de falla 
desconocida. 



Figura 5. Sismo Calderón estación AEPN vs GMPEs 
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ATKINSOM&BO0RE 200B vs SEÑAL SELECCIONADA 






- Med-Geométrica Registro 

Boore &Atki nson2Q0£ 












0 12 3 4 

, m 


Atkmson & Boore 

t 2008 

Mw(S-hB) = 

P 51 

Rjb (0-200 km) = 

27.5 

Vs3J^130-l3ífflms)= 

|~45i:i 

Vref (m/s)= 

P 760 

Tipo de Falla = 

Reverse v 

pga_low= 

P 0.09 

PGAref. (cnVs2H 

100.25 

PGA (cm/s2)= 

34.20 


Figura 6. Sismo Calderón estación LILI vs GMPEs 


ATKIN50N&B00RE JOCB vs SEÑAL SELECCIONADA 



JM 


Figura 7. Sismo Calderón estación PRAM vs GMPEs 


En las figuras 5, 6 y 7, se muestra la comparación de los espectros de 
respuesta, que mejor se ajustan al modelo de predicción de Atkinson&Boore 2008, 
los eventos están dentro de la desviación estándar, denominados “Sismo Calderón 
estación AEPN”, “Sismo Calderón estación LILI” y “Calderón estación PFtAM”. 

La ubicación de las estaciones “AEPN”, “LILI” Y “PRAM”, están dentro del 
Distrito Metropolitano de la ciudad de Quito-Ecuador. 
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6 MODELO DE BINDI ET AL. 2013 (B13) 

Este modelo de predicción del movimiento fuerte (GMPEs), fue desarrollado 
para Europa del Este, derivado de un fuerte banco de datos, especialmente de Italia 
(HACA) y Turquía. (RESORSE DATA BANK), Incluyen 5882 señales, de 1814 
sismos ocurridos en esta zona, registrados entre 1967 hasta el 2011, con rangos de 
magnitudes entre 2.8 y 7.8. 

Los parámetros que calcula este modelo es la aceleración pico del terreno 
PGA, PGV y espectro con la pseudo-aceleración al 5% de amortiguamiento vs el 
periodo, entre 0.002 y 3Seg. 

La forma funcional del modelo es: 

log 10 Y = e 1 + FD(R, M) + FM (M) + FS + F sof (3) 


Y, está dado en cm/s 2 . 

Donde la función de la distancia FD, está dada por la siguiente formula: 

FD(R, M ) = [ Cl +c 2 {M- M ref )\ * log w (Vi? 2 + A 2 /R ref ) - c 3 (Vi? 2 + ft 2 - R ref ) (4) 


La función de la magnitud FM está dada por: 


t6 3 (M-M„) 


para M <M h 
para M > M h 


(5) 


Para FS: 


FS = ylog 10 (Vs 30 /V ref ) (6) 

Donde V re f=800m/s, determinado por regresión, según modelo de Bindi et al. 2013. 


6.1 Características del modelo de Bindi Et Al. 2013. 


El rango de magnitud (Mw) entre 4 a 6.9; distancia a la fuente entre 0- 
200km, estilos de falla validos: Normal, Reverse, Strike-slip y falla desconocida., 
rango del periodo: 0.02-3.0 seg y el tipo de suelo se clasifica de acuerdo al 
Eurocódigo EC8. 

El modelo cuenta con cuatro tipos fuentes sismogénicas de tipo normal con 
593 registros, inversa con 87 registros, strike-slip con 61 registros y falla 
desconocida 28 registros. Además el cálculo de las componentes mediante el 
promedio aritmético de los dos registros seleccionados (mean). 
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GMPEs BINDI 2013 vs REGISTRO ACELEROGRAFICO 




P'ericido(s) 


Bindi13 PlotLog Espectro-Inter, Espectro-Loc. 


Figura 8. Interface programa GMPRs-EPNvl, para el modelo de Bindi et al. 2013. 


En la figura 8, se muestra el espectro de respuesta calculado con el 
programa GMPE-EPN.vl, en el programa, obtenida con la ecuación de atenuación 
propuesta por Bindi et al. (2013). De acuerdo a lo siguientes datos: Mw=6, distancia 
F!¡b=10km, velocidad de onda de corte Vs3o=760m/s y falla tipo strike slip. 

En las figuras 9,10 y 11, se realiza la comparación de los sismos locales 
registrados en Quito y modelo de Bindi et al 2013. Con las señales más cercanas al 
modelo de predicción. 



Figura 9. Sismo Calderón estación AEPN (a) 














































































































218 


Fernando Cañizares, Juan Carlos Singaucho 



Figura 10. Sismo Calderón estación LILI (b) 



Figura 11. Sismo Calderón estación PFiAM 
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7 MODELO DE AKKAR &BOMMER 2013 

Es un modelo actual de movimiento fuerte para la predicción del espectro 
elástico de respuesta en términos de pseudo aceleraciones espectrales, PGA y 
velocidad, usando la base de datos de Italia y Europa, este modelo incluye, la 
amplificación no lineal del suelo en función de la velocidad de onda de corte Vs30, 
con aceleración referencial pico en la roca (Vref). 

El rango de magnitud de momento entre 4 hasta 7, así como, la distancia 
de aplicación cubre los 200km de extensión, el rango de periodos de 0.01 a 4seg. 
Considera diferentes modelos de distancia a la fuente. 

La forma funcional de predicción del movimiento fuerte es: 
lnOO = ln[Y REF (M w ,R, SoF]) + \n[S(y S30 ,PGA REP )] + £a (7) 


7.1 Características del modelo (AB13) 

La Magnitud valida está en el rango de Mw: 4 a 7.6 grados, tipo de fallas: 
normal, inversa, strike slip, distancia Rjb entre 0 - 200km, profundidades menores 
a 30 km. (Depth), para sismos corticales la clasificación del suelo de acuerdo al 
Eurocódigo: Vs30 (150-1200) m/s. 

El modelo incluye el efecto de amplificación no lineal de sitio en función VS30, 
las respuestas obtenidos en términos de PGA, PGV y periodos en el rango de 
T=0.01- 4seg. La definición de la componente horizontal e mediante la media 
geométrica de los registros en las estaciones. 


AKKAR & BOMMER 2013 vs SEÑAL SELECCIONADA 
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Figura 12. Modelo de Akkar&Bommer 2013, GMPE’s EPNvl 
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En la figural 2, se muestra el espectro de respuesta calculado con el 
programa GMPE-EPNvI. Con la ecuación de atenuación propuesta por 
Akkar&Bommer (2013), de acuerdo a los siguientes datos: Magnitud Mw=5, 
distancia Rjb=30km, velocidad de onda de corte Vs3o=800m/s. 


AKKAR& BOMMER2013 vs SEÑAL SELECCIONADA 



0 12 3 4 

_ m _ 

Figura 13. Sismo Calderón estación AEPN. 



Figura 14. Sismo Calderón estación LILI 
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AKKAR& BOMMER 2013 vs SEÑAL SELECCIONADA 



Figura 15. Sismo Calderón estación PFiAM 

En las figuras 13,14y 15. Se muestra los espectros de respuesta entre los 
sismos que mejor se ajustan a la ecuación de atenuación son: Calderón estación 
AEPN, Calderón estación PRAM, Calderón Replica estación AEPN, Calderón 
Replica estación PRAM, corroborando el mismo comportamiento con los modelos 
de Atkinso&Boore 2008 y Bindi et al 2013. 


8 ANÁLISIS DE RESULTADOS ENTRE GMPEs Y SEÑAL 
REGISTRADA 

El análisis de resultados se realiza mediante el índice numérico estadístico 
log-likelihood, y representa la comparación entre dos funciones, el grado de 
cercaníay similitud entre dos funciones f(x) y g(x), en este caso la relación entre el 
espectro de los sismos registrados(fx) y el espectro de la ecuación de atenuación 
propuesta g(x), el cual nos servirá para visualizar el comportamiento de la 
aceleración en los periodos (PGA, 0.1 , 0.5 , 1 y 2 segundos), que finalmente nos 
servirá para recomendar la mejor ecuación de atenuación . 
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Figura 16. índice log-likelihood (LLH) para Atkinson & Boore 2008 


La figura 16, representa el cálculo del índice numérico log-likelihood (LLH) 
en el programa Matlab para el modelo de Atkinson&Boore 2008. La función de 
densidad de probabilidad de color “marrón”, significan los espectros de respuesta 
obtenidos con la ecuación de atención propuesta (GMPEs), en tanto que la otra 
función de densidad de color “azul” representan los espectros de respuesta de las 
señales registradas o seleccionadas, este índice resulto: LLH-PGA= 3.9173, es 
decir, es la diferencia logarítmica entre las dos funciones de densidad de 
probabilidad. Este mismo proceso se sigue para los periodos: 0.1, 0.5, 1 y 2 
segundos respectivamente. 


Los histogramas de color “azul oscuro”, es una ayuda del algoritmo de 
programación para representar las variables de las funciones de densidad de 
probabilidad respectiva en el intervalo de interés. 

Una manera de verificar que el modelo estadístico sea adecuado, es que en 
los periodos bajos (PGA y 0.1 s), este índice debe ser mayor; LLH=3.1050 (los 
modelos GMPEs, no predicen adecuadamente en periodos bajos, en relación al 
valor real registrado), en tanto que el índice LLH debe ser menor; LLH=2.2702, para 
valores de periodo alto (2 segundos). 
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Figura 17. índice log-likelihood (LLH) Bindi 2013 


El mismo análisis de la figura 16, se lo realiza para las figuras 17 y 18, el 
cual representa el cálculo del índice LLH, para los modelos de Bindi et al (2013) y 
Akkar&bommer 2013. Para los periodos de interés (PGA, 0.1, 0.5, 1 y 2 segundos) 


3 Figure 1 — □ X 

File Ed'rt View Insert Tools Deslctop Window Help 

Da a ¿ fe \ % f; © ^ y. -1 a □ §hj o 



z z z 



z z 


Figura 18. índice log-likelihood (LLH), para Akkar&Bommer 2013 































































224 


Fernando Cañizares, Juan Carlos Singaucho 



Figura 19. índice log-likelihood de los GMPEs parasismos corticales 

La figura 19, se obtiene graficando los pedidos de interés (PGA, 0.1, 0.2, 
0.5, 1 y 2 segundos, mostrados en la columna 1 de la Tabla 3, con los respectivos 
valores calculados de LLH (Columnas 2,3 y 4 de la Tabla 3), para los diferentes 
modelos de atenuación y se generan una líneas de tendencia para poderlos 
comparar gráficamente. 


Tabla 3. Valores log-likelihood de las GMPEs para sismos corticales. 


PERIODO 

ÍNDICE LIKELIHOOD (LLH) PARAGMPE s Y REGISTROS CORTICALES 

Atkinson&Boore 2008 

Bindi 2013 

Akkar&Bommer 2013 

0.001 

3.9173 

3.1052 

3.9627 

0.1 

3.8651 

3.2199 

3.7475 

0.5 

3.3376 

2.5781 

2.9951 

1 

2.8108 

2.2557 

2.4318 

2 

2.5649 

2.2702 

2.4682 


El menor índice log-likelihood (LLH), mostrado en la figura 19, y en la Tabla 
3, es el modelo Bindi et al. (2013). Por tant, este modelo es el que mejor se ajusta 
entre los datos reales y los calculados, sin embargo, los modelos de Atkinson 
&Boore (2008) y Akkar& Bommer (2013), también muestran una buena tendencia. 


9 CONCLUSIONES 

Es importante el cálculo de la longitud y forma de la falla, mediante las 
ecuaciones de Wells y Coppersmith (1994), (Rjb). Este dato influye en la forma de 
las respuestas esperadas, debido a que las fuentes sismogénicas corticales están 
cerca de la ciudad de Quito (5-25Km). 

De acuerdo al modelo basado en el análisis estadístico log-likelihood (LLH), 
mostrado en la figura 23, el menor índice de ajuste entre los valores reales y 
calculados, es la ecuación de atenuación de Bindi et al. (2013), para sismos de 
origen cortical. Por lo tanto, se podrá utilizar este modelo para el análisis del riesgo 
sísmico para la ciudad de Quito. 
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Se recomienda que para una futura actualización de la norma Ecuatoriana 
de Construcción (NEC 15), o para estudios de riesgo sísmico de Quito o de ciudades 
con características geológicas similares, se pueda utilizar los modelos de 
atenuación de Akkar y Bommer 2013 y Bindi 2013. Los mismos que incorporan 
nuevos parámetros en las ecuaciones como: la amplificación lineal y no lineal del 
suelo, además se cuenta con la validación de los eventos reales registrados en la 
red RENAC, del Instituto Geofísico de la Escuela Politécnica Nacional (2015). 

Los modelos seleccionados son aplicables para sismos corticales entre 
rangos de magnitud “bajos” (Mw=4-5), como los registrados en la ciudad de Quito. 
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